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Abstract 

In 1872 G. Darboux denned a family of curves on surfaces of M 3 
which are preserved by the action of the Mobius group and share many 
properties with geodesies. Here we characterize these curves under the 
view point of Lorentz geometry and prove some general properties and 
make them explicit them on simple surfaces, retrieving results of Pell 
(1900) and Santalo (1941). 

Introduction 

Our interest here is to understand a family of curves on a surface called 
Darboux curves. Almost as in the case of geodesies, through every point and 
direction which is not of principal curvature, passes a unique Darboux curve. 
References about these curves are [Dal] . [Ri] . [Co] . [Pej . [En] . [Salj . |Sa2j . 

To understand the dynamics of these curves we will give a natural family 
of spheres along a Darboux curve in the set of spheres of 1R 3 or S 3 . Then 
we will study the angle drift of Darboux curves with respect to the foliations 
of the surface by principal curvature lines and by curves making a constant 
angle with the principal curvature foliations. 

We will also study the Darboux curves on some special surfaces such 
as the envelopes of very particular one-parameter family of spheres, and on 
quadrics. 
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This paper is organized as follows: In section [T] is described the space of 
spheres of S 3 and the correspondence with a quadric Acl 5 with the induced 
Lorentz metric of signature 1. In section [2] a comparative study of quadrics 
and Dupin cyclides is studied considering a foliation by constant angle with 
the principal lines. In section [3] is the introduced the space of osculating 
spheres to a surface and the Darboux curves are characterized geometrically. 
In section 2] is performed a geometric study of Darboux curves on Dupin 
cyclides. In section [5] is obtained the differential equation of Darboux curves 
in a principal chart and this is applied in subsequent sections. 

In section [6] is defined a natural plane field in the space of spheres and the 
study of its integrability is carried out. In section [7] a study of Darboux curves 
near regular ridge points leading to the zigzag and beak to beak behaviors 
is carried out. In [8] is considered the study of Darboux curves in general 
cylinders, cones and surfaces of revolution, viewed as canal surfaces. In 
section [9] is carried out the global study of Darboux curves on quadrics. 

1 The set of spheres in S 




Figure 1: S 3 ^ and the correspondence between points of A 4 and spheres. 

The Lorentz quadratic form C on M. 5 and the associated Lorentz bilinear 
formC(-, •), are defined by C(x , ■ ■ ■ , x 4 ) = — x"l + (x\ + - ■ ■+x 2 4 ) and C(u, v) = 
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-UqVq + {UiVi H h U4VA). 

The Euclidean space IR 5 equipped with this pseudo-inner product C is 
called the Lorentz space and denoted by 1L 5 . 

The isotropy cone Li = {v G M 5 | £(t>) = 0} of £ is called the light cone. 
Its non-zero vectors are also called light-like vectors. The light cone divides 
the set of vectors v G 1L 5 , v {L = 0} in two classes: 

A vector v in M. 5 is called space-like if £(i>) > and time-like if < 0. 
A straight line is called space-like (or time-like) if it contains a space-like 
(or respectively, time-like) vector. 

The points at infinity of the light cone in the upper half space {xq > 
0} form a 3-dimensional sphere. Let it be denoted by S^,. Since it can 
be considered as the set of lines through the origin in the light cone, it 
is identified with the intersection of the upper half light cone and the 
hyperplane {x = 1}, which is given by S 3 = {(xi, • • • ,x 4 ) | x\ + - ■ -+x\ — 1 = 



To each point a G A 4 = {v G M 5 | C(v) = 1} corresponds a sphere 
£ = a 1 - fl or E = er 1 fi (see Figure [[]). Instead of finding the points 
of S 3 "at infinity" , we can also consider the section of the lightcone by a 
space-like affine hyperplane H z tangent to the upper sheet of the hyperboloid 
7i = {C = — 1} at a point z. This intersection Light fl H z inherits from the 
Lorentz metric a metric of constant curvature 1 (see |H-Jj . |La-Wa] . and 
Figure Wj. 



Figure 2: A tangent space to M 5 cuts the light cone at a unit sphere 

Notice that the intersection of A 4 with a space-like plane P containing 
the origin is a circle 7 C A 4 of radius one in P. The points of this circle corre- 



0}. 



A 
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Figure 3: Spherical and Euclidean models in the Minkowski space 1L 5 (up). 
The geodesic curvature k g , picture in the affine hyperplane TC (down). 



spond to the spheres of a pencil with base circle. The arc-length of a segment 
contained in 7 is equal to the angle between the spheres corresponding to 
the extremities of the arc. 

It is convenient to have a formula giving the point a G A 4 in terms of the 
Riemannian geometry of the corresponding sphere E C S 3 C Light and a 
point m on it. For that we need to know also the unit vector n tangent to S 3 
and normal to £ at m and the geodesic curvature of £, that is the geodesic 
curvature k g of any geodesic circle on S. 

Proposition 1. The point a G A 4 corresponding to the sphere £ C S 3 C 
Light is given by 

a = k g m + n. (1) 

Remark: A similar proposition can be stated for spheres in the Euclidean 
space E 3 seen as a section of the light cone by an affine hyperplane parallel 
to an hyperplane tangent to the light cone. 

The proof of Proposition [T] can be found in [H-Jj and [La-Ohj . The idea 
of the proof is shown on Figure [3J Let 7i be the affine hyperplane such that 
S 3 = Light fl H, let P be the hyperplane such that £ = S 3 PI P; the vertex 
of the cone, contained in 7i, tangent to S 3 along E is a point of the line P 1 - 
which contains the point a G A 4 . From Proposition [[J we see that the points 
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Figure 4: Tangent spheres 



in A 4 corresponding to a pencil of spheres tangent to a surface Mat a point 
m form two parallel light-rays (one for each choice of normal vector n) . Let 
us now chose the normal vector n, and consider the spheres £& associated to 
the point o> = km + n. All the spheres but for at most two, have either 




Saddle contact 



L'Oiitact with an 
osculating sphere 



Center contact 



Figure 5: Possible contacts of a sphere and a surface 



a center contact or a saddle contact with M (see Figure [5]). 

The exceptional spheres correspond to the principal curvatures of M at 
m, they are called osculating spheres. 

When k ^ [fci,/^] the intersection of £& and M near the origin reduces 
to a point, the origin. 

When k E (ki, k-i) the intersection of £& and M near the origin consists 
of two curves intersecting transversely at m. 

When k = k\ or k 2 the intersection of £& and M near the origin is a 
singular curve, in general of cuspidal type, at m. 

In fact one can prove the following 
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Proposition 2. Let k = kicos 2 a + k2sin 2 a. Then the angle of the tangents 
at m to Sfc fl M with the principal direction corresponding to k\ is ±a. 

In general, the intersection of an osculating sphere with M admits a 
cuspidal point at m, the tangent to the cusp is then the principal direction 
associated to the curvature fcj. In any case, when k G [ki, k 2 ] goes to k iy the 
two tangent directions at m to fl M converge to the principal direction 
associated to the curvature ki. From that we can see that lines of curvature 
as osculating spheres are conformally defined. 

2 Foliations making a constant angle with re- 
spect to Principal Foliations 

Consider a surface M with principal foliations V\ and P 2 and umbilic set U. 
The triple V = (Pi, Vi,W) will be referred as the em principal configuration 
of the surface. 

Definition 3. For each angle a G (— 7r/2,7r/2) we can consider the foliations 
JF+ and T~ such the leaves of this foliation are the curves making a constant 
angle ±a with the leaves of the principal foliation V\ . We will write T a = 



In other words, the normal curvature of a leaf of T a is precisely k n {a) = 



Figure 6: Foliation of a torus of revolution by Villarceau circles 



Dupin cyclides are very special: they are surfaces which are in two dif- 
ferent ways envelopes of one-parameter families of spheres (see |Da3j ). This 



k\ cos 2 a + k 2 sin 2 a. 



2.1 Foliations T a on Dupin cyclides 
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implies that the corresponding curves are circles or hyperbolas in A 4 , inter- 
section of A 4 with an affine plane (see [La-Waj ). 

There are three types of Dupin cyclides. One can chose a nice representant 
of each class: 

- A) The boundary of a tubular neighbourhood of a geodesic of S 3 

- B) A cylinder of revolution in R 3 

- C) A cone of revolution in R 3 . 

Then, in cases A) and B) the foliations and T~ are totally geodesic 
foliations. 

In the case B), four foliations are foliation by circles: the two foliation by 
characteristic circles JF and J r 7T /2, and two others: the foliations by Villarceau 
circles. 

In the case C), we have to develop the cone on a plane -this is a local 
isometry out of the origin- to see that a foliation of the plane by curves 
making a constant angle with rays is a foliation by logarithmic spirals. The 
picture on the cone is obtained by rolling the plane foliation back on the 
cone. 



2.2 Foliations T a on quadrics 

2 2 2 

Proposition 4. Consider an ellipsoid E aj b ]C = {{x, y,z) : — + 4- + — = 1} 
with a > b > c > 0. Then E a 6 iC have four umbilic points located in the 
plane of symmetry orthogonal to middle axis; they are of the type D X) i.e., a 
singularity of index 1/2 and having one separatrix for the principal curvature 
lines. For all a the singularities of T a are the four umbilic points and this 
configuration is topologically equivalent to the principal configuration V of 
the ellipsoid near the umbilic points which are of type D\ 

Proof. Consider the parametrization of the ellipsoid in a neighborhood of the 
umbilic point p = (u , 0, v ) = (yj^jf, 0, yj £ ^ ). 



a(u,v) = p +uE 1 +vE 2 +-[Vac(u 2 +v 2 )+^—^ -(u 3 +uv 2 )+h.o.t]E 3 +h.o.t. 

2 b 6 

Here {Ei, E 2 , E 3 }, E 2 = (0,1,0), is a positive orthonormal base and the 
ellipsoid is oriented by E 3 = —VH(p )/\VH(p ). 

In a neighborhood of the umbilic point (0, 0) the differential equation of 
the foliation T a in the chart (u, v) is given by: 

A(u, v)dv 2 + B(u, v)dudv + C(u, v)du 2 = 0, where 
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A(u, v) = — u — cos 2a \/ u 2 + v 2 + R 3 (u, v) + A 2 (u, v) 
B(u,v) =2v + B 2 {u,v) 

C(u, v) =u — cos2a\Ai 2 + v 2 + R 3 (u,v) + C 2 (u, v) 

Here A 2 = 0(r 2 ), C 2 = 0(r 2 ) and R 3 (u,v) = 0(r 3 ), r = y/u 2 + w 2 . 

The above implicit equation has, real separatrices with limit direction 
given by ±2a and the behavior of the integral curves near is the same of 
an umbilic point of type D\. 

In fact, consider the blowing-up u = rcosa, v = rsina. 

The differential equation of T a in the new variables is given by: 

( cos 2a — cosa + rR 1 (r, a))dr 2 + r(2 sina + rR 2 (r, a))drda 
+r 2 (cos 2a + cos a + rR 3 (r, a))da 2 = 0. 

The two singular points are given by (2a, 0) and (—2a, 0). Direct analysis 
shows that both singular points are hyperbolic saddles of the adapted vector 
fields to the implicit equation near these singularities. The blowing-down of 
the saddle separatrices are the umbilic separatrices of JF+ and T~. See Fig. 

El 




Figure 7: Resolution of the foliations JF+ and T a . 

Therefore it follows that the pair of foliations and T~ near an umbilic 
point of the ellipsoid with three distinct axes is topologically equivalent to 
the configuration of principal lines near a Darbouxian umbilic point D±. □ 

Remark: The study of principal lines was first considered by Monge, see 
|Molj and |Mo2] . Near umbilic points the behavior of principal lines on real 
analytic surfaces was established by Darboux, [Da2] . See also [Guj . |G-S] . 
|Ga-S3j and references therein. 

Proposition 5. Consider an ellipsoid E a] b iC with a > b > c > 0. On the el- 
lipse Yi xz C E a ,ft,c> containing the four umbilic points, Pi, (i = ,4 ) 
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counterclockwise oriented, denote by Si(ct) = 2 J c sina LjwvJ rftt ( resp. 

52(a) = 2 J fc a cosa[-j^j]dv ) a distance between the adjacent umbilic points p\ 

and pi ( resp. p\ and P2 )■ Define p(a) = f4^y- 

T/ien if p <ER\Q ( resp. p G Q) all the leaves of ' T a are recurrent ( resp. 
all, with the exception of the umbilic separatrices, are closed). See Fig. [21 




Figure 8: Foliations IZi of the ellipsoid E a fe c 



Proof. The ellipsoid E a ft jC belongs to the triple orthogonal system of surfaces 
defined by the one parameter family of quadrics, -^j + + = 1 with 
a > b > c > 0, see also |Sp| and [St]. The following parametrization a(u, v) = 
(x(u, v),y(u, v), z(u, v)) of E aM , where 



/ a(u - a)(v - a) b(u-b)(v-b) c{u-c){v-c) \ 
a(U ^ )= n (b-a)(c-a) (b-a)(b-c)^]l (c-a)(c-b) J' (2) 

defines the principal coordinates (u, w ) on E a 6 iC , with u G (6, a) and i> G (c, 6). 
The first fundamental form of E a ,b,c is given by: 

I = ds 2 = Edu 2 + Gcfe 2 = v - — { du 2 + v — -{ dv 2 3 

4ff(u) AH(v) 

The second fundamental form with respect to the normal N = — (a u A 
a v )/\a u A a„| is given by 



(v-u) labc.n (u-v) abc, 2 



4ff(u) V uv 4H(v) V w (4) 

H(x) =(x — a)(x — — c) 
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Therefore the principal curvatures are given by: 



e 1 abc g 1 abc 
h = — = -\ , k 2 = — = -\ . 

h, u V uv G v V uv 



The four umbilic points are given by: 



The differential equation of the foliation T a in the principal chart (u, v) 
is given by 

H(u)v cos 2 a dv 2 + H{v)u sin 2 a du 2 = -v^ 

2 i 2 '2 ^ i 2 r\ 

cos a dv + sin a ——^du = 0. 



ifM H (u) 



Define <io"i = sin a ^_^^ du and g?c72 = cos a j^yjdu 

Therefore the differential equation of T a is equivalent to da\ — do\ = 
in the rectangle [0, si(a)] x [0, S2(a)J. 

On the ellipse E = {(x,y, z)\{— + — = 1, y = 0} define the distance 
between the umbilic points pi = (x , 0, 2 ) and p 4 = (x , 0, — 2 ) by si(a) = 
2 J c fe sin a[j^0^}du and that between the umbilic points p\ = (xo, 0, zq) and 

= {-Xo, 0, z ) is given by s 2 (a) = 2 J^[cosa^}dv. 

The ellipse E is the union of four umbilic points and four principal um- 
bilical separatrices for the principal foliations. So E\{j>i,p 2 ,P3,P4} is a 
transversal section of the foliation T a - 

Therefore near the umbilic point p\ the foliation, say JF*, with umbilic 
separatrix contained in the region {y > 0} define a the return map a + : E — ► 
E which is an isometry, reverting the orientation, with a + (pi) = p\. This 
follows because in the principal chart (u, v) this return map is defined by 
er + : {u = b} — > {v = b} which satisfies the differential equation ^ = — 1. 
By analytic continuation it results that a + is a isometry reverting orientation 
with two fixed points {pi, P3}. The geometric reflection cr_, defined in the 
region y < have the two umbilic {p 2 , Pa} as fixed points. So the Poincare 
return map 7Ti : E — > E (composition of two isometries a + and cr_) is a 
rotation with rotation number given by S2(a)/si(a). 

Analogously for J 72 with the Poincare return map given by tt 2 = r + o r_ 
where r + and r_ are two isometries having respectively {p 2) Pi\ and {pi,P3} 
as fixed points. □ 
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Remark: The special case a = ir/4 was studied in [Ga-Slj . A more gen- 
eral framework of implicit differential equations, unifying various families of 
geometric curves was studied in [Ga-S2j . See also |Ga-S3j . 

Proposition 6. In any surface, free of umbilic points, the leaves of and 
T~ are Darboux curves if and only if the surface is conformal to a Dupin 
cyclide. 

Proof. From the differential equation of Darboux curves, see equation (fTUl) 
in Section [5], it follows that: 

nn , , . da 1 dki 3 1 dk 2 . 3 

6{ki — K2) smacosa —— = —=—— cos a H -=— — sin a. 

ds y/E du VG dv 

So all Darboux curves are leaves of if and only if ^(w, v) = ^ (u, v) = 0. 
By Proposition [25] in Section [8] this is exactly the condition that characterizes 
the Dupin cyclides. □ 



3 Spheres tangent to a surface along a curve 

Above each point m e M which is not an umbilic, the spheres tangent to 
the surface M having a saddle contact with M form an interval of boundary 
the two osculating spheres at m. Let us first define the 3-dimensional subset 
V(M) C A 4 as the union of the spheres having a saddle contact with the 
surface M and the osculating spheres. 

Therefore V(M) is an interval fiber-bundle TrV(m) —>■ M over M; the 
boundary of V(M) is the surface O of spheres osculating M (this surface 
has two connected components when the surface M has no umbilical point). 
V(M) is a 3-manifold with boundary a surface O made of spheres osculating 
M (this surface has two connected components when the surface M has no 
umbilical point). V(M) inherits from A 4 a semi-Riemannian metric. At each 
point, the kernel direction is the direction of the light ray through the point. 

Let us consider now a curve C C M 2 . 

The restriction of V(M) to C form a two-dimensional surface in A 4 which 
is a light-ray interval bundle V(C) over C out of the umbilical points of M 
which may belong to C. From V(M), we get on V(C) an induced semi- 
Riemannian metric. 

In this text we will use for derivatives with respect to parametrization 
of curves contained in R 3 or S 3 , (often the parameter is an arc length), and • 
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for derivatives with respect to a parametrization of a curve in A 4 (often the 
parameter is an arc length, but now for the metric induced from the Lorentz 
"metric"). 

One section of this bundle has the property that at every point m where 
the curve is not tangent to a principal direction, the point o~c>(m) correspond 
to a sphere E(m) such that one branch of the intersection Ep/(m) fl M is 
tangent to C at m. If at m the curve C is tangent to a principal direction 
we take the corresponding osculating sphere as Ec'(m). 

We call this section C ansec{C) of V(C) the canonical section, and the 
corresponding family of spheres the canonical family along C; the envelope 
CanCan(C) of the spheres of this family is called canonical canal correspond- 
ing to C C M. 

Proposition 7. Consider the notation above. Then the following holds: 

i) The section Cansec(C) satisfies C(k g ) > 0, and moreover k g e T m V(M). 

ii) The section Cansec(C) is a geodesic in V(C). It is of minimal length 
among sections ofV(C). Moreover it is the unique section on V(C) which 
is of critical length. 

Hi) One of the cuspidal edges of the envelope CanCaniC) is C. 

Proof. The condition defining the sphere E^(m) implies that the order of 
contact of C and Ec(m) is at least 2, one more that the order of contact of 
Ec/(m) H M) and C, which is at least one. To verify this property, notice 
that, in terms of the arc-length of a branch of Ec(m) fl M), or equivalently 
of the arc- length s on C, the angle of E(m) with M along Ecy(m) fl M) is 
of the order of s, if not smaller. The distance to C is of the order of s 2 , if 
not smaller. Therefore the order of the distance of a point of C to Ec(m) 
is of order s 3 if not smaller. This means that C and the sphere E^/ (m) have 
contact of order at least 2 at m. 

Definition 8. A curve on M such that at each point its osculating sphere is 
tangent to M will be called Darboux curves. 

Recall a lemma from [L-Sj 

Lemma 9. A curve T(t) = span(7(t)) has contact of order > k with a sphere 
E corresponding to a iff 

a±span( 7 (t), 7 (t),..., 7 ( fc )(t)) 
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Proof. The sphere £ is the zero level of the function f(x) = (x\o~). Then the 
contact of T(t) and £ has the order of the zero of (/ o j)(t) = (-f(t)\a). □ 

Therefore, the sphere Cansec(m) is orthogonal to m, m and m. This 
ends the proof of Proposition UJ □ 



Differentiating the relation (er|m) = and using the relation (o~\fh) = 
we get (a\m) = 0. Differentiating the relation (a\m) = 0, and using the 
relation (<r|m) = 0, we get the relation (a\m) = 0. 

We can use as parameter the arc-length of the space-like curve a = 
Cansec. Its geodesic curvature vector, orthogonal to rh is therefore orthog- 
onal to m © K. a, proving that {a} is a geodesic in V(C). 

We notice that a is tangent to T m V(M), as T m (V(M)) = T m (Cight) = 
m ± . 

Now we will use the Darboux frame T, Ni,n,m of the curve C C M C 
S 3 C Cight, where Ni is the unit vector tangent to M normal to C compatible 
with the orientation of M. Using again the formula a = k n m + n, where 
k n is the normal curvature in the direction tangent to C at m, we see that, 
when a is the section Cansec(C) 

\a'\ = \k' n m + k n m' + n' — \k' n m — r g Ni\ = \r g \ (5) 

where the geodesic torsion r g is defined by the following formula 

((V T n) | N x ) = -r g . (6) 

Observe that our formula ([S]) gives an interpretation of the geodesic torsion 
of a curve Cc Mas the rotation speed of the canonical family of spheres 
tangent to M along C. 

In order to compute the geodesic curvature vector of the section Cansec(C), 
we need to use its parametrization by arc- length in A 4 . Then a = a — = 

-^{k' n m — TgNi). Differentiating once more, we get 



r 1 
a = — 9 -{k' n m- TgNi) + -^\k" n m + k' n T - t'N\ - r g (-k g T + r g n)}, 



T 3 Vn 9 u ' 1 T 2< 

9 9 



notice that k g is the geodesic curvature of C C M, while k g is the geodesic 
curvature vector of the curve Cansec C A 4 . 
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Simplifying, we get 



a = cj)(s)m + — [(k' n + T g kg)T - r 2 g n] 
T g 

As kg — a + a we get 

kg = ip{s)m + ^(k' n + Tgkg)T (7) 

As the vector k g is orthogonal to m and to a = —(k' n m — r g Ni), we see 
that the curve Cansec is a geodesic on V(C). 

Otherwise, when a = km + n, we see that spheres tangent to a surface 
along a curve form a space-like curve in A 4 ; explicitly we get 

\a\ = |-[(A;- k n )'m + (k - k n )T + Tg N x \\, (8) 

giving, as m, T and N are mutually orthogonal, a proof of the fact that the 
section Cansec(C) has minimal length among sections, as the three vectors 
m,T,N are mutually orthogonal in L 5 . Formula [8] shows also that no other 
section of V(C) is of critical length. 

Remark: The characteristic circle of the envelope is the intersection of 
span(a, a) ± and the sphere S 3 C 1R 4 ; as the vector T is orthogonal to m, n 
and Ni, and therefore to a and a' (and also a), the characteristic circle is 
tangent at m to C. 

Recall that a drill ( |Thoj . |L-Sj ) is a curve in the space of spheres the 
geodesic curvature of which is light-like at each point. Generically, points of 
a drill are osculating spheres to the curve C G S 3 defined by the geodesic 

acceleration vector k g of the drill. 

We see that if we can find drills in V(m) we find geodesies. In fact we 
find that way almost all of them. 

Formula [7] implies, as m is the normal direction to V(M) at a point 
a = km + n, the following theorem. 

Theorem 10. The curve C is a Darboux curve if and only if the section 
Cansec(C) C V(M) is a geodesic in V(M). This happens if and only if the 
curve C C M satisfies the equation: 

K + Tgk g = 0. (9) 
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Remark: The light rays of V(M) are also geodesies. Segments of geodesies 
of V(M) which are not tangent to light rays define an arc of curve C C M, 

and therefore, when k g is light-like, C is a piece of Darboux curve on M. 

Also when C is a Darboux curve the only cuspidal edge of the envelope 
CanCan(C) is C. See jTho] . 

Definition 11. We denote by E mj £ the sphere tangent to the surface M at 
the point m such than one branch of the intersection E • DM is tangent to 
the direction I (canonical sphere associated to the direction £). We will also 
use the notation E m >v when the direction £ is generated by a non-zero vector 
v. 

Given a curve C C M such that the tangent vector to C at c(t) is con- 
tained in I , E,^ • is, among the spheres tangent to M at m, the one which 
have the best contact at m with the curve C . 

Remark: r g ds is the differential of the rotation of the sphere E c , t s along 
the curve C. 

Proof: The sphere E c ^ has a tangent movement which is a rotation of 
"axis" the characteristic circle of the family, which is tangent to C . It is 
tangent to the tangent plane. We see that changing the sphere tangent to 
the surface along C changes the "pitch" term but not the "roll" term. □ 

Proposition 12. Let C be a curve contained in the surface M C S 3 . Let 
7c be the curve in A 4 obtained considering at each point m £ C the sphere 
S m £( m ) ; where l(m) is the direction tangent at m to C . Suppose that the 
curve C is nowhere tangent to a principal direction of curvature. Then the 
curve 7c is space-like with at every point a space-like geodesic acceleration. 
Moreover the curve C is one fold of the singular locus of the canal surface 
defined by jc- 

Proof: The condition defining the canonical sphere E,^ • ^ guarantees that 
the characteristic circle CC(t) is tangent to C, which is therefore the singular 
curve of the canal. □ 

Proposition 13. The condition of the equation in the theorem [TO) is equiv- 
alent to the fact that one branch of the intersection of the sphere E,^ . has 
at the point the same geodesic curvature as the curve C . 

Proof: As the sphere E c ^ • is osculating the curve C when it is a Darboux 
curve, the contact order with the curve should be at least 3. □ 
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When the direction i defined by c(t) is not a principal direction, there is 
a unique sphere tangent to M which has contact with C of best order, the 
other having contact of order one with C. There is only one such sphere, in 
the pencil containing the osculating circle to C which is tangent to M. This 
last sphere should therefore be E,^ - , t -. 

Remark: The previous considerations imply that the geodesic curvature of 

—k f 

the branch of S a fl M tangent to the Darboux curve is equal to: k g = -^ JL . 

4 Darboux curves in cyclides: a geometric 
viewpoint 

We can describe V(M) when the surface M is a regular Dupin cyclide (see 
Figure [9]). It is the wedge of the two circles formed by the osculating spheres 
of the regular cyclide. 



Figure 9: Darboux curves in V(M) when M is a regular Dupin cyclide 

Consider a surface M with principal foliations V\ and V2 and umbilic 
set U. The triple V = (Vi,V2,M) will be referred as the em principal 
configuration of the surface. 

Definition 14. For each angle a G (— 7r/2,7r/2) we can consider the folia- 
tions JF+ and T~ such the leaves of this foliation are the curves making a 
constant angle ±a with the leaves of the principal foliation V\ . 

In other words, the normal curvature of a leaf of T a is precisely k n (a) = 
k\ cos 2 a + k,2 sin 2 a. 

Let us, for later use, define the surfaces M a C V(M) as the set of spheres 
tangent to M at a point m having curvature k a = k n (a) = ki cos 2 a + 




Osculating Spheres (fei) 



mean spheres 



Osculating Spheres (fc 2 ) 
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&2 sin 2 a. Each surface M a is foliated by the lifts of the curves of T a \ we call 
this new foliation T a . These later foliation form a foliation of V(M) that we 
call T. 

Proposition 15. The Darboux curves in Dupin cyclides are the leaves of the 
foliations and T~ . 

The proof of Proposition [T51 is given in subsections 14.0.11 and 14.0.21 
4.0.1 Cylinders in M 3 

This is the easiest case to visualize. Cylinders are models for cyclides having 
exactly one singular point. 

The definition of a Darboux curve requires that the osculating sphere to 
the helix 7i drawn on the revolution cylinder be tangent to it. This comes 
from the fact that the helix is invariant by the rotation R of angle tt and 
axis equal to the principal curvature vector of the helix which is a vector 
orthogonal to the cylinder and going along its axis of revolution. As the 
osculating sphere should also be invariant by this rotation its diameter in 
contained in the line normal to the cylinder which is the axis of the rotation. 
Therefore the sphere is tangent to the cylinder. 

4.0.2 Regular cyclides in S 3 

Each regular Dupin cyclide M is, for a suitable metric of constant curvature 
1, the tubular neighborhood of a geodesic Q\ of S 3 (see |La-Waj ). Seeing S 3 as 
the unit sphere of an euclidian space E 4 of dimension 4, it is the intersection 
of a 2-plane Pi with S 3 . Let P<i be the plane orthogonal to Pi in E 4 . Then it 
is also a tubular neighborhood of the geodesic Q 2 = P2 H S 3 . We can define 
(for this metric) the symmetries with respect to the two spheres containing 
respectively m and Q\ and m and Qi- Let 1Z be the composition of these 
symmetries. Then 1Z preserves the cyclide but also the "helices" on the 
cyclide. It should therefore also preserve the osculating sphere to the helix, 
which, as it is not normal to the cyclide, has to be tangent to it. The "helix" 
is therefore a Darboux curve. There are enough "helices" to be sure we got 
all the Darboux curves. 

The last family of Dupin cyclides is formed of conformal images of cones 
of revolution (see [La-Wa] ). On can deal with cones of revolution as we 
did with regular cyclides, using a sphere orthogonal to the axis and m (it 
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belongs to the pencil whose limit points are the singular points) and a sphere 
containing the axis and m. 

We will, using the general dynamical properties of Darboux curves, prove 
(see Proposition [6]) that a surface is a Dupin cyclide if and only if its Darboux 
curves are the leaves of the foliation T of V(M). 



5 Differential Equation of Darboux Curves in 
a Principal Chart 

Consider a local principal chart (u,v) in a surface M C M 3 . The first and 
second fundamental forms are denoted by 

I = Edu 2 + Gdv 2 , 11 = 1 = edu 2 + gdv 2 

and the principal curvatures are k\ = e/E and k 2 = g/G. 

Proposition 16. Let (u, v) be a principal chart. Let c be a curve parametrized 
by arc length s making an angle a(s) with the principal direction (1,0). For 
a Darboux line c the following differential equation is verified 



da 1 dki , 1 dk 2 o 

3(ki — k 2 ) sm a cos a — — =—=— — cos a H ;=-^— sin a 

ds ^E du VG dv 

dki n dk 2 . q 
cos a + — — sm a. 



(10) 



ds i ds 

Here {v! ,v') = (^3,^3). 



Proof. Consider a principal chart (u, v) such that v = cte are the leaves of 
the principal foliation V\. 

Let c(s) = (u(s),v(s)) be a regular curve parametrized by arc length s. 
So we can write c'(s) = (u',v') = ("^fiT^f)) defining a direction a with 
respect to principal foliation Vi, horizontal foliation. 

We have the following classical relations: 



k n {a) =k\ cos a + k 2 sin a 

kq -f k\ cos a + k 2 sin a 
ds y y 

T g =(k 2 — ki) cos a sin a. 
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Here k g = {k g ) v=const and k 2 g = (k g ) u=const are the geodesic curvatures of 
the coordinates curves. 
Therefore, 

dk n 1 dk\ o 1 dk\ 2 1 ^2 2 

= ^=— — cos cH — = — — cos a sin cH — =— — cosasm a 

as y/E ou ov y/E ou 

1 8^2 . s n/i 7 \ .da 
H — sin a + 2(/c2 — «a) cos asm a— 
VG <9v (is 



The differential equation of Darboux curves is given by k' n + k g r g = and 
so it follows that: 

[kl(k 2 — fci) H — 7=^7^] cos 2 a sin a + [A; 2 (/c2 — &i) H — ^"7r~] cos a sin 2 a 
y/G ov y J E ou 



da 1 dki □ 1 <%2 q 
+ 3{ko — k\ ) cosasma- — I = cos a -\ -=— — sm a = 

In any orthogonal chart (F = 0) we have that G u = 2G\/~Ek 2 g = (k g ) u=cte ) 
and E v = —2E\fGk x g = {k g ) v=cte . Also the Codazzi equations in a principal 
chart are given by: 

dh E v dk 2 G u 

^7 = 2E ik2 - kl) > ^v- = 2G {h - h) - 

See Struik's book [53 pages 113 and 120]. 

Therefore, 



^-*>> + ^t= s^-^S* -*>=»■ 

This ends the proof. □ 

Remark: A curve c(s) has contact of third order with the associated oscu- 
lating sphere, tangent to the surface, when 

(c',c')[2(N',c") + (N",c')} - 3(c',N')(c',c") = 0. 

This equation can be used to obtain the differential equation of Darboux 
curves in any chart (u, v). See |Salj . 
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6 A plane-field on V(M) 

The two tangents to the two Darboux orbits through the point (m, a) G 
V(M) define a plane in T m ^ a V(M). The ensemble of these planes define a 
plane-field V . 

Proposition 17. The 'plane-field V is integrable if and only if 

(6)02 = -^0102, (6)01 = ^0102- 

O 

Here 6 are the conformal vector fields and 9i are the principal conformal 
curvatures. 

Proof. Consider in the unitary tangent bundle the suspension of Darboux 
curves. 

These curves are defined by the following vector field 



cos a d sin a d dk% / du cos 2 a dk 2 / dv sin 2 a d 

1 \[E du s/G dv 3VE(ki - k 2 ) sin a 3y/E(ki - k 2 ) cos a da' 

Consider the involution <p(u, v, a) = (u, v, —a) and the induced vector 
field D 2 = <^(£>i). 
So it follows that: 

cos a d sin a d dki/du cos 2 a dk 2 /dv sin 2 a d 

2 \[E du s/G dv ?,VE{kx-k 2 ) sin a 3y/E(h - k 2 ) cos a 1 da' 

Consider the plane field, Darboux plane field, defined by {Di,V 2 }. 
In the sequence it will be obtained the conditions of integrability of this 
plane field. 

In order to simplify the calculations the following changes will be devel- 
oped. 

First, consider the new pair of vector fields T>\ — T>\ + T> 2 and T> 2 = 
T>\ — T> 2 and obtain: 

~ 2 cos a d r 2dk-i/du cos 2 a, d 

1 y/E du 3v A E(A;i - k 2 ) sin a 1 da 
~ 2 sin a d . 1dk 2 ldv sin 2 a, d 

2 VG dv [ 3VE(h-k 2 ) cos a 1 da' 
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Next consider: 

- 2 d . 2dk 1 /du cos a, 9 

1 x/Edw 3V^(fci - fc 2 ) sin a J 9a 

- 2d, 2dk 2 /dv sin a, 9 
25 = l r li l 

2 x/G^ 3-/^(fci - fe) cos a da 

Consider the unitary vector fields Xi, the conformal vector fields £j and 
the principal conformal curvatures 0j. 



X 1 

6 



l 9 



2X l 



ki — k 2 ' 
, 4(X 1 )fc 1 



X 2 = 

6 = 
o 2 



dki/du 



Observing that — H= 
base defined by: 



1 d 
V&dv 
2X 2 
h - k 2 
4(X 2 )fc 2 

(h - k 2 

8k2/dv 



.2 ■ 



X i (k i )=Dk i (X i ). 



(Xi)ki and H^- = (X 2 )k 2 we obtain a new 



v c 2 =6 



1 cos a 9 

6 1 sin a da 

1 sin a 9 

6 2 cos a 9a 



Recall that: 



[fX,gY] = fg[X,Y] + (X.g)fY - (Y.f)gX 

[£i,6] ^ 

So it follows that: 



--02 fl - 2 61 ! ^2- 



[£ £ l _|_ [£ 1 ^ sin « 9 1 ^ cos a 9 
' ' 6 cos a da 6 sin a 9a ' 



1 cos a 1 sin a 

+ c^ 1 " — Tr^® 2 TT\ 

b sin a oa b cos a oa 

_ - - i 9lf2 + i( fl )0^ £ - i<6*!55H #■ 

2 2 b cos a oa b sm a oa 



1 



9 



W-iflifc + kfOfc— -(6)0: 
2 2 b cos a 



cos a 



1 



18 sin a cos a 9a 
102 , ^ 



sin a 3 sin a cos a da 
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Now consider the determinant 



1 






1 



sin a 
sin a 



A 



1 



[(£l)02 



sin a 



(6)0 



cos a 



1 0i0 2 
3 sin a cos a 





A 



6 



cos a 



sin a 



Evaluating this determinant it is obtained: 



-3[(£i)02 + cos 2a + [3(6)^2 - 3(6)0! + 0i0 2 ] 



36 sin a cos a 
So the integrability conditions are given by: 



(6)^2 + (6)01 = 0, 

This is equivalent to: (6)02 



3(6)02-3(6)0i + 0i02 = o. 

-\e x Q 2 and (6)01 = ±0i0 2 . 



□ 



7 Darboux lines near a Ridge Point 

Consider a surface M and a principal chart (it, v) such that the horizontal 
foliation V\ is that associated to the principal curvature h\. 

Definition 18. A non umbilic point p Q = (u , Vq) is called a ridge point of 
the principal foliation V\ if 7^7 G°o) = 0. 

The ridges for the principal foliation V2 are characterized by ^f-(po) = 0. 

Definition 19. A ridge point relative to V\ is called zigzag, respectively beak 
to beak, when ai(p ) = ^r{po)/{h{po) ~ hi'Po)) < 0, respectively, ai(p ) > 
0. 

A ridge point p Q relative to V2 is called zigzag, respectively beak to beak, 
when a 2 (p ) = -gjr(Po)/ (h(Po) - hipo)) < 0, respectively cr 2 (Po) > 0. 

The ridge points are associated to inflections of the principal curvature 
lines, to the singularities of the focal set of the surface and also with the 
singularities of the boundary of V(M), the space of spheres tangent to M. 
See |Poj for an introduction to ridges and also [Guj in his study of geometric 
optic and applications in construction of eye lens. 

A practical way to see the type of a ridge point is given by the following 
proposition. 
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Proposition 20. Consider a surface of class C r , r > 4 parametrized by the 
graph (u, v, h(u, v)) where 

, . k\ 9 9 CI q d 9 b n Co 

v ' ; 2 2 6 2 2 6 

A , B o C 99 D n E A 

H it H it f H w w H w H u + /i.o.t 

24 6 4 6 24 



T/ien (0, 0) is a ridge point for V\ when a = 0. Also o~i = [ ^J 3 ^ 1 + t^- 



2f/ 2 



(fci-fc 2 ) 2 J - 

Corresponding to the foliation P 2 ^ e point (0, 0) zs a ridge point when 
c = 0. ^^(0) = ^ + ^]. 

Proof. Straightforward calculations shows that the principal curvatures in a 
neighborhood of (0, 0) are given by: 



fci(it, v) =ki + au + dv + \{A - 3kf + 2d )u 2 + (B - 2 — )uv 

2 K\ — K2 R2 — k\ 

1 2b 2 

+ (C - hk 2 2 - - -)v 2 + h.o.t. 

2 K2 — k\ 

k 2 (u, v) =k 2 + bu + cv + \{C - k\k 2 + ^ )m 2 + (-D + 2 M )uv 

1 26 2 

+-(E - 3k* + - -)v 2 + h.o.t. 

2 k 2 — ki 

(11) 

So the result follows. □ 

Proposition 21. Let p be a ridge point of M corresponding to principal 
foliation V\ such that o~i(po) ^ 0. Then the ridge set R containing p is locally 
a regular curve transversal to V\ and the boundary of V(M) corresponding 
to a = has a cuspidal edge along ii~ l (R). See Fig. [JJJ Analogously for the 
ridges associated to the principal foliation P 2 - 

Proof. In the parametrization given in Proposition [2D] it follows that the at 
Po = (0,0) the principal direction corresponding to V\ is t\ = (1,0) and 
the ridge set is parametrized, according to equation (fTTj) . by: (u(v),v) = 

Next consider a principal chart (u,v). The set of spheres V(M) C 
A 4 is parametrized by o~(u,v,a) = k n (a)m(u,v) + N(u,v) with k n (a) = 
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Figure 10: Ridges and Singularities of the boundary of V(M) 

k\{u, v) cos 2 a + k 2 (u, v) sin 2 a, £(m, m) — and 7V U = —kim u , N v = —him. 
See equation (JTj). We have that 



So -D0 has rank 3 for a G (0, |). 

The boundary of V(M) is parametrized by a = and a = ir/2 and so 
0i (it, u) = ki(u, v)m(u, v) + N(u, v) and 02(1*, v) = k 2 {u^ v)m(u, v) + N(u, v). 
The map 01 has rank 1 at the ridges and so we have the structure of cuspidal 
edges on the boundary of V(M). □ 

Proposition 22. Let po be a non umbilic and disjoint from the ridge set. 
The Darboux curves tangent to the principal lines are given by cuspidal curves 
= (lfr( fc i - h)* + ■ ■ ■ . IK** " k ^ 3 + ■■■) and r 2 (t) = (ff (fc 2 - 
ki)t 3 + ■ ■ ■ , \^{k 2 — ki)t 2 + ■••). The behavior of all the Darboux curves 
passing through p is as shown in the Fig. [711 

Proof. Consider the vector field X defined by the differential equation: 




dh 2 dk 2 . 2 1 , /, 
o~ u = \-^r~ cos a + — sm aim + (k ri 
ou ou 



dk, 



a v cos a + — — sin aim + (k n 

ov ov 

o~ a —\{k 2 — k\) cos a sin aim 
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Figure 11: Darboux curves through a non ridge point 



u =—= cos a[3(ki — k 2 ) sin a cos a] 
v E 



v' =—= sin a[3(/ci — k 2 ) sin a cos a] 



. 1 o 1 9A; 2 o 
a = —=— —cos cH =— — sin a. 

The projections of the integral curves of X in the coordinates (u, v) are 
precisely the Darboux curves. 

For any initial condition (0,«o), with a ^ kir/2, the integral curves 
of X are transversal to the axis a and so has a regular projection. For 
«o = n7r/2 and a 2 ^ 0, the integral curves of X are tangent to the axis 
a and the projections are of cuspidal type. For a = nn direct calculations 
gives: (u(t), v(t)) = (fff (A* - fc 2 )t 2 + ■ • • , (-lff^i - ^ 2 )t 3 + •■•)• Now 
observe that the projection of the integral curves passing through (0, 0, 0) 
and (0, 0, 7r) are both tangent to semi axis of u. □ 

Theorem 23. Let R be an arc of ridge points transversal to the corresponding 
principal foliation, i.e., suppose that o~i(p) 7^ for every p G R. Then there 
are two types of behavior for the Darboux curves near the ridge set, the first 
is the zigzag and the other is the beak to beak. 



Proof. Consider the vector field X as in the proof of proposition [22j 

We will consider only the ridge set corresponding to V\. For the other 
principal foliation the analysis is similar. The ridge set is defined by the 
equations v) — and ^ (u, v) = 0, one for the corresponding principal 

foliation. 
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Figure 12: Darboux curves near regular curve of ridges: zigzag and beak to 
beak. 



The singularities of X are defined by (U (v), v), 0) and (u, V(u), |), where 
fH^),^) = and ff K V>)) = 0. 

To simplify the notation suppose a singular point (0, 0, 0) of the ridge 
transversal to the principal foliation V\ and E(0) = G(0) = 1. 

It follows that: 

/ 3(h-k 2 y 
DX(0) =00 
V& o 
The eigenvalues of DX(0) are: 



1 d 2 ki 1 d 2 ki 

Al = °' A2 = 7! V ~d^ /{kl ~ hl A2 = "7! V ~d^ /{kl ~ k2) - 

By invariant manifold theory, when A2A3 = — fe) = — |°"i(0) < 

0, the singular set of X (ridge set) is normally hyperbolic and there are 
stable and unstable surfaces, normally hyperbolic along the singular set. This 
implies that there is a lamination (continuous fibration) along the ridge set 
and the fibers are the Darboux curves. Also the prolonged Darboux curves 
are o class C 1 along the ridge set. 

So the Darboux curves are as shown in Fig. [121 center and right. That 
is, there are Darboux curve crossing the ridge, tangent to the principal lines, 
and the prolonged Darboux curves are C 1 along the ridge set. 

In the case when ^(0) = ^-/(h -k 2 )<0 the non zeros eigenvalues of 
DX(0) are purely complex and so the singular set is not normally hyperbolic. 
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In this case we are in the hypothesis of Roussarie Theorem, [Rol Theorem 
20, page 59], so there is a local first integral in a neighborhood of the ridge 
set. The level sets of this first integral are cylinders and the integral curves 
(helices) in each cylinder when projected in the surface M has a cuspidal 
point exactly when helix cross the section a = 0. This produces the zigzag. 

There are no Darboux curves tangent to the principal direction e\ along 
the ridge set in this case. □ 

8 Darboux curves on general cylinders, cones 
and surfaces of revolution 

Darboux curves on general cones where already studied by Santalo ( |Sa2j ). 
In a similar way, one can study Darboux curves on cylinders and surfaces of 
revolution. This is not a coincidence. The three type of surfaces are canal 
surfaces corresponding to a curve 7 C A 4 which is also contained in a 3- 
dimensional subset of lA Depending on the subspace, this intersection is 
either a copy of A 2 , a unit sphere S 3 or a 2-dimensional cylinder (see [Dal] , 
[M-Nj . |Ba-La-Wa] l The latter condit ion defines conformal images of general 
cones, general cylinders and surfaces of revolution. 

These surfaces can be obtained imposing conformally invariant local con- 
ditions. 

Recall that, assuming that S is a surface which is umbilic free, that is, 
that the principal curvatures k\{x) and k2(x) of S are different at any point 
x of S. Let X\ and X 2 be unit vector fields tangent to the curvature lines 
corresponding to, respectively, ki and k 2 . Throughout the paper, we assume 
that k\ > k 2 . Put \i — (k± — k 2 )/2. Since more than 100 years, it is known 
([Tr]. see also [CSWj ) that the vector fields & = Xi/ fi and the coefficients 
(i = 1,2) in 

[6, 6] = ~(0 2 fi + 0i6) 

are invariant under arbitrary (orientation preserving) conformal transforma- 
tion of R 3 . (In fact, they are invariant under arbitrary conformal change 
of the Riemannian metric on the ambient space). Elementary calculation 
involving Codazzi equations shows that 

e l = \-X 1 {k 1 ) and 6 2 = \-X 2 {k 2 ). 
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The quantities 9i (i — 1,2) are called conformal principal curvatures of S. 

Canal surface are characterized locally, |Dalj . |H-Jj . |M-Nj . |Ba-La-Wa] . 
by the following propositions. 

Proposition 24. A surface S is (a piece of) of a canal if and only if one of 
its conformal principal curvatures, say 9 2 , is equal to zero. 

Proposition 25. Imposing moreover that the other conformal curvature, say 
6\, is constant along characteristic circles characterizes the surface as one of 
the three families above (cone, cylinder and surfaces of revolution) . 

Proposition 26. Let M be a surface and (u, v) be a principal chart such that 
9i(u,v) = 6\{u) and 9 2 (u,v) = 0. Let A(u) = exp[J k ^l k , 2 du] and a G (0, ir) 
be an angle. Then the function J{u, a) = A{u) cos 3 a is a first integral of the 
Darboux curves. Moreover in the region A c = ir(M c ) = {(u, v) : u e M c }, 
M c = l 7~ 1 (c), the Darboux curves are defined by the implicit differential 
equation 

c^Gdv 2 - E(A 2/Z - c 2 ' z )du 2 = 0. 
Proof. The differential equation ( flQj) is simplified to 

. cos a . sin a . 1 k[ cos 2 a 

U = V = Ct — ; — 

4E ' VG ' 3v^ h - k 2 sin a ' 

1 1 k 1 

So it follows that 5 s = ;1 — l_£2L^ which is a equation of separable vari- 

du A fci— K2 sin a ^ 1 

ables. Direct integration leads to the first integral J as stated. 

To obtain the implicit differential equation solve the equation J{u, v ) = c 
in function of cos a and observe that 4^ = ^Ssmo^ q 

du i/G cos a 

Proposition 27. Let M be a surface and (u, v) be a principal chart such 
that 9i(u, v) = 9\{u) and 9 2 (u, v) =0. Then the plane-field V is integrable. 

Proof. Direct from the characterization of integrability of V established in 
Proposition \T7\ □ 

8.1 Darboux lines on surfaces of Revolution as Canal 
Surface 

Consider an one parameter family of spheres of radius r(u) with center at 
(0,0, u). 



28 



The envelope of this family is a canal surface and can be parametrized 

by: 



H(u, v) = r(u) cos/9(w)(cosf , sinv, 0) + (0, 0, u — r(u) sin P(u)), 

where cos f3{u) = — r'(u) 2 , sin/5 = r', \r'(u)\ < 1, (3 E (-7r/2, tt/2). 
The normal unitary to the surface is: 

N = (— cos (3{u) cos v, — cos sinv, sin f3{u)). 
The coefficients of the first and second fundamental forms of H are given 



by: 



E(u, v) = (1 f f ^ ) \ F{% v) = 0, G(«, v) = r 2 (l - r' 2 ) 
1 — r 



e(«,v) = ^ , /(«,v)=0, = r(l - r' ). 

1 — r 

The principal curvatures are given by: 

r" 1 
h(u, v) = — -2 fc 2 («, u) = -. 

It will be assumed that k 2 > &i and the surface is free of umbilic points. 
The ridge set is defined by the equation 

Bin, v) = ^ = r"'(l - r' 2 ) + 3rV' 2 = 0. 
ou 

Proposition 28. The Darboux lines on surfaces of revolution, free of umbilic 
points, can be integrated by quadratures. The function 



T(u, a) = r cos (3{kx — k 2 ) cos 3 a = h(u)(ki — k 2 ) cos 3 a, 



is a first integral of the differential equation of Darboux lines. Here h is the 
distance of the point of the surface to the axis of revolution. 

Moreover, if R'(u) < the ridge is zigzag. If R'(u) > the ridge is beak 
to beak. See Fig. \TT 
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Figure 13: Darboux lines near regular ridges on surfaces of revolution 



Proof. In the principal chart (u, v) the differential equation of Darboux lines 
is given by 

. 1 .1 . 1 k[ sin 2 a 

u = — = cosa, v = — sma, a 



VE VG 3 VE h - k 2 1 cos a 

rs that 
Now observe that 



So it follows that 3^da = -j^du. 

cos a k\ — ki 



J ^« = / f^f *< + / ^* = Hh -*.)+/ ^tr/u 

2 ——du = \ 1 7y]du = lnr(l — r' ) 

h-k 2 J r 1-r' 21 

Therefore it follows that 

Z(w, v, a) =r(l — r' 2 ) l l 2 {ki — k 2 ) cos 3 a 

=r cos/3(/ci — k 2 ) cos 3 a = h{u){ki — k 2 ) cos 3 a 

is the first integral. The analysis of behavior near ridges follows from Theo- 
rem E3 □ 

Proposition 29. The Darboux lines on a cone, free of umbilic points, can 
be integrated by quadratures. The function 

Z(w, a) = k g {u) cos 3 a 

is a first integral of the differential equation of Darboux lines. Here k g is the 
geodesic curvature of the intersection of the cone with the unitary sphere. 
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Figure 14: Darboux lines near regular ridges on a general cone 



Moreover, if k' g /k g < the ridge is zigzag. If k' g /k g > the ridge is beak 
to beak. See Fig. \T^\ 

Proof. The cone can be parametrized by X(u,v) = vy(u) where |7| = 1 and 
|7'| = 1 is a spherical curve. We have that ki(u, v) = k g {u) and k 2 (u,v) = 0, 
since 7" = —7 + k g ^ A 7' and N(u, v) — 7 A 7'. 

The Darboux curves are given by: (sinct/ cosa)<ia = ^(k' g /k g )du. So 
it follows that X(u,a) = k g (u) cos 3 a is the first integral. The analysis of 
behavior near ridges follows from Theorem [23j □ 

Proposition 30. The ridges on a surface of revolution are characterized by 
the singularities of the curve 0\ C A 2 C A 4 defined by the osculating circles 
of envelope of plane curves with center (0,u) and radius r(s). 

Moreover the ridges correspond to "vertices" of the curve 7 C A 2 C A 4 
defining the surface of revolution as an envelope. 

Proof. The envelope of plane curves with center (0,u) and radius r(u) is a 
plane curve c(u) with curvature given by k{u) = — 1 _ r ,,a_ rT „ • 

According to equation flTJ the curve Oi(u) = k{u)m{u) + n{u) with 
C(m,m) = 0, C(n, n) = 1 and n'(u) = —k{u)m'{u). So it follows that 
0[ = k'{u)m{u) is a light vector ant is zero precisely at the points where 
k'\u) = 0. 

Let us remark that if a point of 7(5) 6 7 is a vertex then the radius p(s) 
of the osculating circle Os) is critical. Let us denote by P(s) the osculating 
plane to 7 at 7(5), and by x(s) G P(s) the center of the osculating plane, 
which is also the (unique) critical point of the restriction of £ to P(s). The 
Dupin necklace formed by "the other osculating spheres" along the charac- 
teristic circle C(s) C S(s), where the sphere S(s) correspond to the point 
7(5) e A 4 is the envelope of the spheres of the circle Os) and of the sphere 
of the associated circle 0*s). This circle is contained in the affine plane 
Q(s), orthogonal to span(0, P(s) and which contains the point y(s) which 
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belongs to the line spanned by x(s) and satisfies £(x(s),y(s)) = 1. The 
square Lorentz norm C(y(s )) is critical when C(x(s )) and p(s ) are. Let 
a(s) be a curve of osculating spheres of the other family along a line principal 
curvature. It is a light-like curve. Notice that, as Q(s) is orthogonal to y(s), 
C(y(s),a(s) = C(y(s)). Derivating £(y(s),a(s)) with respect to s we get, 
for a critical s , g)| s=so + £( J, <r(s))| a=ao = 0. As y(s) belongs to 

the 3-dimensional space A = span(Os),0), which is independant of s, ^| 
belongs also to A. All the affine spaces Q(s) are orthogonl to A, therefore 
£(%ct(s))\ s=so = C(%y(s))\ s=so = 0. Therefore C(y(s), f )\ s=so = 0. This 
is possible if and only if ^\ S=SQ = 0, as a Lorentz scalar product of non-zero 
time-like and light-like vector cannot be zero. The sphere is therefore a 
ridge sphere. This reasonning is valid for all the curves in A 4 formed of os- 
culating spheres of the other family along a line of principal curvature. The 
whole characteristic circle C (so) is therefore a ridge. □ 

Remark: Similarly in a cylinder a(u, v) = c{u) + v z , where c is a plane 
curve with curvature k the function T(u, a) = k{u) cos 3 a is a first integral 
of Darboux curves. 

Remark: The geodesies on surfaces of revolution has a first integral given 
by J~(u,a) = h(u)sina. The parallel u = Uq is a geodesic if and only if 
h'(uo) = 0, and it is a hyperbolic geodesic when h"(u ) > 0. In a cone a first 
integral for the geodesies is given by J{v, a) = cosa/v. 

9 Darboux curves on quadrics 

The Darboux curves in the ellipsoid were considered in [Pej by Pell. Here we 
complete and hopefully simplify his work. 

The quadrics Q a ,b,c belongs to the triple orthogonal system of surfaces 
defined by the one parameter family of quadrics, + + = 1 with 
a > b > c > 0, see also [Sp] and [St] . 

Consider the principal chart (u, v) and the parametrization of Q a ,b,c given 
by equation (j2J). 

For the ellipsoid u G (6, a), v G (c, b) or u G (c, b), v G (b, a). 

For the hyperboloid of one sheet u G (b,a), v < c or u < c, v G (b, a). 

For the hyperboloid of two sheets u G (c,b), v < c or u < c, v G (c,b). 

The first fundamental form of Q a ,b,c is given by equation ([3]) and the 
second is given by equation (jll). 
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Therefore the principal curvatures are given by: 



e 1 abc g 1 abc 

h = — = -\ , k 2 = — = -\ . 

E u V uv G v V uv 

Also consider the functions 

ri = fl/(3( fcl - fa)) = i^JL_, r2 = ^ /3 ((*, _ * 2 » = I. « 



9m 2u(m-jj)' <9t> 2 d(u - u) 



a(a-b) n — I — / c(6-c) 



The four umbilic points are (±xo, 0, ±2 ) = (±y °_ c , 0, ± 

Proposition 31. TTie differential equation of Darboux curves on a quadric 
Q a ,b,c is given by: 

u =—= cos 2 a sin a, u' = — == cos a sin 2 a 



E VG (12) 

a = — cos =sin a 

y/E VG 

Proof. In a principal chart the differential equation of Darboux curves is 
given by equation (jTUJ). Define tana = ^gj so that Eu' 2 + GV 2 = 1. To 
obtain a regular extension of the differential equation (fTOj) to a = mir/2 and 
consider it as a vector field in the variables (u, v, a) multiply the resultant 
equation by the factor cos a sin a and the result follows. □ 

Proposition 32. The function 



ti i a ri \ cos a sin a v G(u,v) , . 

l(u, v , u , v ) — I(u, v, a) — 1 , tana = — W — (13) 

u v u' V E(u, v) 

is a first integral of equation ( fI2i) . 
Proof. We have that I(u, v,a) 



We will show that 4-(I(u(s),v(s), a(s)) = along a solution (u(s),v(s), a(s)). 
We have that 

1 ?/??' 1 ?/' 

1 ??7/' 1 ?J 

Straightforward calculation leads to 4-(I(u(s),v(s),u'(s),v'(s))) = 0. 

□ 
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Proposition 33. The Darboux curves on a quadric Q a ,b,c ore the real integral 
curves of the implicit differential equation: 

(v - X)H(u)v' 2 - (u - X)H{v)u' 2 = 0. 

The normal curvature in a Darboux direction T> defined by I(u, v, dv/du) = 

1/A is given by k n (p, V) = 

This differential equation is equivalent to 



e(u,v)du 2 + g(u,v)dv 2 1 abc 1 1/4 1/4 

Proof. From the equation I — 1/ A it follows that 

,dv. 2 (u — X)vE(u,v) (u — X)H(v) 



du (v — A) uG(u, v) (v — X)H(u) 
Here dv/du is a Darboux direction T>. 

As k n = EXctdvfdu) 2 ^ f° nows from the equation above and from equations 

and Q that k n {p,V) = \^ = \{abc) 1/A K}/\ K = hk 2 . 

For the reciprocal part consider the implicit differential equation 



e(u,v)du 2 + g(u,v)dv 2 1 abc 
E(u,v)du 2 + G(u,v)dv 2 X V uv 

Using equations Q and (j4j) it follows that this equation is equivalent to 
the following. 



-—-du 2 —dv 2 = (u-X)H(v)du 2 -(v-X)H(u)dv 2 = 0. 

H(u) H{v) 

The restriction in the values of A is in order to consider only real solutions 
of the implicit differential equation obtained. □ 

Proposition 34. The ridge set of the quadric Q a ,b,c is the intersection of the 

quadric with the coordinates planes. Moreover: 

a) For the ellipsoid with 0<c<b<ait follows that: 

i) The ellipse E xy = {z = 0} PI ^ a ^ c , respectively E yz = {x = 0} PI ~E, a ,b,c, 
is a ridge corresponding to k2, respectively to k\, and is zigzag. 
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ii) The ellipse E xz = {y = 0} H E a ^ c containing the four umbilic points 
(±x , 0, ±z ) is the union of ridges of ki and k 2 . For \x\ > x the ridge 
correspond to k±. The ellipse E xz is beak to beak in both cases. 

b ) For the hyperboloid of one sheet with c<0<b<ait follows that: 

i) The hyperbole H yz is a ridge corresponding to k\ and is beak to beak. 

ii) The hyperbole H xz is a ridge corresponding to k\ and is zigzag. 
Hi) The ellipse E xy is a ridge corresponding k 2 and it is zigzag. 

c) For the hyperboloid of two sheets with c<b<0<ait follows that: 

i) The hyperbole H xz is a ridge corresponding to k\ and is zigzag. 

ii) The hyperbole H xy containing the four umbilic points (±xo,±2/o?0) is 
the union of ridges of k\ and k 2 . For \x\ > xq the ridge correspond to 
k\ and all segments of hyperbolas are beak to beak. 

Proof, a) Ellipsoid: By symmetry is clear that the points of intersection 
of the coordinates planes with the ellipsoid are ridges points. The principal 
curvatures have no critical points along the corresponding principal curvature 
line in the complement of these three ellipses. In a principal chart (u, v) we 
have that = -f k x ^ and = -f k 2 ^ 0. 

It will be sufficient to check the condition of zigzag or beak to beak only 
in a point of a connected component of the ridge set. 

Consider the point po = (~\A^ 0? 0)- The ellipsoid is parametrized by: 

x(y, z) = -Va- + yft- + - + — + — + — + h.o.t] 

Therefore ki(p ) = y/a/b, k 2 (p Q ) = y/a/c, A = 3y/a/b 2 (A is the 
coefficient of y A ) and ((A - 3kf)(k 2 - h) = -3a(a - b)(b - c)/(6 4 c) < 0. 
Therefore the ellipse E xz is beak to beak. 

Also, let E = 3y/a/c 2 (coefficient of z A ). So, (E — 3k 2 )(k\ — k 2 ) = 3a(a — 
c)(b — c)/(6c 4 ) > 0. Therefore, by Theorem [23] the ellipse E xy is zigzag. 

Now consider the point go = (0, — Vb, 0). The ellipsoid is parametrized 

by: 

•JC^ 1 rJC^ z^ 1 z^ 1 

y(x, z) = -Vb + Vb[— + — + — + - — + — + h.o.t.] 
v ' l 2a 2c 8a 2 Aac 8c 2 J 
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Now, with A coefficient of x A and E coefficient of z A it follows that (A — 
3k\){k 2 - h) = 3b(a -b)(a- c)/(a A c) > and (E - 3k^)(k 1 - k 2 ) = 3b{b - 
c)(a- c)/(ac 4 ) > 0. 

So both ellipses E yz and E xy are both zigzag, one for the corresponding 
principal curvature. 

b) Hyperboloid of one sheet: the ridges are given by the intersection of the 
hyperboloid with the coordinates planes. 

Consider the point go = (0, — yb, 0). The hyperboloid is parametrized by: 



,-y»2 ~2 /-y>4 , yi 2 ~2 —2 

y(x, z) = -Vb + Vb[— + — + — + - — + — + h.o.t. 

2a 2c 8a 2 Aac 8c 2 



Therefore h(q ) = Vb/a > k 2 (q ) = Vb/c < 0, A = 3y/b/a 2 (A is 
the coefficient of x A ) and E = 3\fbjc 2 coefficient of z A it follows that (A — 
3fc?)(fci - k 2 ) = -3b(a - b)(a - c)/(a 4 c) > and (E - 3Jfc|)(Jfc 2 - h) = 
—3b(b — c)(a — c)/(ac 4 ) < 0. The hyperbole H yz is beak to beak and the 
ellipse E xy is zigzag. 

Next consider the point po = (—y/a, 0,0). The hyperboloid is parametrized 

by: 

( \ r< r\y 2 j_ 2:2 j_ y A j_ ^ 2 j_ ^ 2 j_ a + i 

Therefore fci(po) = v^/^ > k 2 (p ) = y/a/c < 0, A = 3y/a/b 2 
(A is the coefficient of y A ) E = 3yfajc 2 (coefficient of z A ) it follows that 
((A - 3kf)(h - k 2 ) = 3a(a - 6)(6 - c)/{b A c) < and (E - 3k%)(k 2 - h) = 
—3a(a — c){b — c) / \bc A ) < 0. Therefore the hyperbole H xz and the ellipse E xy 
are both zigzag. 

c) Hyperboloid of two sheets: the ridges are the intersection of the hyper- 
boloid with the coordinates planes. 

Consider the point po = (y/a, 0, 0). One leaf of the hyperboloid is parametrized 

by: 

Therefore h(p ) = -y/a/b > k 2 (p ) = -y/a/c > 0, A = -3y/a/b 2 (A 
is the coefficient of y A ) and ((A — 3kf)(ki — k 2 ) = 3a(a — b) (a — c)/{b A c) < 0. 
Therefore the hyperbole H xz is zigzag. 

Also, let E = -3yfajc 2 (coefficient of z A ). So, (E - 3kf){k 2 - fci) = 
— 3a(a — c)(b — c)/{bc A ) > 0. Therefore the hyperbole H xy for \x\ < x is beak 
to beak. 
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For |x| > Xq on the hyperbole H xy a similar analysis shows that it is beak 



Proposition 35. Consider the ellipsoid E ai b jC with a > b > c > 0. 

i) For c < A < b the Darboux curves are and contained in cylindrical 
region c < v < A and the behavior is as in Fig. [731 upper left. 

ii) For A = b the Darboux curves are the circular sections of the ellipsoid. 
These circles are contained in planes parallels to the tangent plane to 
E a ,b. c at the umbilic points. These circles are tangent along the ellipse 
E y and through each umbilic point pass only one Darboux curve. See 
Fig. [T5J, bottom right. 

Hi) For b < A < a the Darboux curves are bounded contained in the two 
cylindrical region A < u < a and the behavior is as shown in the Fig. 
[75l upper right. 



to beak. 



□ 







Figure 15: Darboux curves on the ellipsoid. 



Proof. First case: c < A < b. 

The differential equation of Darboux curves is given by: 




and b < u < a. 



Define do~\ = (u — A)/ H{u)du and do 2 = \J{v — A) /H{v)dv. 
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Therefore the differential equation is equivalent to do~\ — do~\ = , with 
(01,02) G [0,Li] x [0,L 2 ] ( Li = £ da x < oo, L 2 = j c do 2 < oo ). 
In the ellipsoid this analysis implies the following. 

The cylindrical region C\ = a ([b, a] x [c, A]) is foliated by the integral 
curves of an implicit differential equation having cusp singularities in dC\. 
We observe that this region is free of umbilic point and is bounded by prin- 
cipal curvature lines, in coordinates defined by v = c and v = A. 

The case b < A < a the analysis is similar. Now the differential equation 
of Darboux curves are defined in the region [A, a] x [c, b] and we have a 
cylindrical region C\ = a ([A, a] x [c, b}). This 

For A = b the differential equation can be simplified in the following. 

(u — a)(u — c)dv 2 — (v — a)(v — c)du 2 = 0. 

This equation is well defined in the rectangle [c, a] x [c, a]. 

Define d&\ = \j \J (u — a)(u — c)du and da 2 = l/y/(y — a) (y — c)dv. So 
the equation is equivalent da\ — da\ = with (01,02) G [0,L] x [0, L] (L = 
j c dai). So in this rectangle all solutions are straight lines. The images of 
this family of curves on the ellipsoid are its circular sections. In fact we 
know that the ellipsoid has circular sections parallel to the tangent planes at 
umbilic points. As the circles are always Darboux lines it follows that the 
solutions of the differential equation is the family of circular sections. So we 
have two families of circles having tangency along the ellipse E y . □ 

Proposition 36. Consider an ellipsoid ^ a ,b,c with three axes a > b > c> 
and suppose b < A < a. Let L\ := yj E(u, b)du and L 2 := y/ G(b, Xjdu 
and define p = j^. Consider the Poincare map it : S — ► X associated to 
the foliation of Darboux curves defined by the implicit differential equation 
1= 1/A. 

Then if p G M \ Q ( resp. p G Q) all orbits are recurrent ( resp. periodic) 
on the cylinder region v < A}. See Figure [75| bottom left. 

Proof. The differential equation of Darboux curves is given by: 

(V ~ A) ,2 (U ~ A) /2 

— , v T———U =0, c<v<X<u<a. 

H(v) H(u) 



Define d<j\ = J du and da 2 = J %^ dv. By integration, this leads 

to the chart (01,02), in a rectangle [0, L±] x [0,L 2 ] in which the differential 
equation of Darboux is given by 

do\ — da\ = 0. 
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The proof ends with the analysis of the rotation number of the above equa- 
tion. See similar analysis in Propositions and [2SJ □ 

Proposition 37. Consider a connect component of a hyperboloid of two 
sheets M a ^ c with a > > b > c. 

i) For A < c the Darboux curves are non bounded and contained in the 
non bounded region v < A and the behavior is as in the Fig. [Tb\ left. 

ii) For A = c the Darboux curves are the circular sections of the hyper- 
boloid. See Fig. EH center. 

Hi) For c < A < b the Darboux curves are non bounded and contained in 
the cylindrical region A < u < b and the local behavior is as shown in 
Fig. EH right. 




Proof. The analysis developed in the case of the ellipsoid also works here. 
See proof of Proposition [3SJ □ 

Proposition 38. Consider an hyperboloid of one sheet HI aj b jC with a > b > 
> a Let A e (-oo, c) U (b, oo). 

i) For A < c the Darboux curves are bounded and contained in the cylin- 
drical region A < v < c and the behavior is as in Fig. 11, upper left. 

ii) For b < A < a the Darboux curves are unbounded and contained in 
the cylindrical region b < u < A [outside the hyperbola E x ) and the 
behavior is as in Fig. VT% upper right. 



39 



in) For A = c and A = a the Darboux curves are straight lines of the 
hyperboloid. For A = b the solutions are not real. See Fig. [77[ bottom 
left. 

iv) For a < A all Darboux curves are regular ( helices) and goes to oo in 
both directions. See Fig. [77| bottom right. 




Figure 17: Darboux curves on a hyperboloid of one sheet. 



Proof. Similar to the proof of Proposition [35j □ 

Remark: The global behavior of geodesies in quadrics, in particular in the 
ellipsoid, was studied in [Ga-S3j . 
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